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Abstract 

The elastic 7rN scattering amplitude in the isospin limit is calculated in 
the framework of heavy baryon chiral perturbation theory, up to the third 
order. Threshold parameters like scattering lengths, volumes, effective 
ranges, etc. are compared with data. All relevant low energy constants 
are fixed from the available pion-nucleon data. A clear improvement in the 
description of data is observed, when going from the first two orders in the 
chiral expansion to the third one. The importance of even higher orders is 
suggested by the results. 
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1 Introduction 



Constraints of chiral symmetry on pion-nucleon interactions were investigated in 
the sixties in terms of current algebra, and the prediction of the S-wave 7rN scat- 
tering lengths Jl| was one of the most important results within this approach. In 
the eighties, the systematic method — called chiral perturbation theory (CHPT) 
- of calculating corrections to the current algebra results was invented [0], and 
it was applied to 7rN scattering up to order 0(p 3 ) ||. 

However, CHPT with nucleons is not as systematic as CHPT for light mesons, 
since the nucleon mass spoils one of the main virtues of CHPT — one-to-one cor- 
respondence between loop expansion and expansion in external momenta. This 
correspondence is valid for massless particles 0, so it works for the Goldstone 
bosons of the spontaneously broken chiral symmetry of QCD (pions, kaons and 
rj), which are massless in the chiral limit of massless u, d and s quarks. Nonzero 
masses of light quarks, leading to nonzero masses of Goldstone bosons, are treated 
as a perturbation. The simultaneous expansion in external momenta and light 
quark masses is called chiral expansion. Nucleons, on the other hand, are massive 
even in the chiral limit and their masses cannot be treated as small perturbations. 
As a consequence, there is no more direct correspondence between loop and chi- 
ral expansions, and diagrams with arbitrary number of loops can contribute to a 
given chiral order. 

This drawback was eliminated in the clever reformulation of CHPT with 
baryons - called heavy baryon chiral perturbation theory (HBCHPT) - 
which, so to say, shifts the nucleon mass from the propagator to vertices of an 
effective Lagrangian and thus restores the loop-chiral correspondence. HBCHPT 
was applied to different processes up to order 0(p 3 ) || (and even 0(p 4 ) [j7|]|§), 
among others to the forward threshold 7rN scattering, from where the correc- 
tions to the S-wave 7rN scattering lengths up to order 0(p 3 ) were calculated 0. 
Recently also higher 7rN partial waves were calculated up to this order JTDj . 

In these works, some of the low energy constants (LECs) of the 0(p 3 ) ttN 
Lagrangian were not determined from the fit to the 7rN data, instead they were 
estimated from the principle of resonance saturation, which had been shown to 
work very well in the mesonic sector fill , but in the baryonic sector it is just a 



working hypothesis. 

The purpose of this paper is to calculate the full 7rN scattering amplitude in 
the isospin limit in HBCHPT to order 0(p 3 ) and to fix the LECs contributing 
to this amplitude from the available experimental information. The paper is 
organized as follows. In section 2 we summarize the previous results, needed for 
the calculation of the amplitude, in section 3 the amplitude is calculated and in 
section 4 it is confronted with data. Conclusions are summarized in section 5, 
and some technical points and/or lengthy formulae are left to appendices. 
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2 HBCHPT of the pion-nucleon system 
2.1 Effective Lagrangian 

Our calculation of the elastic 7rN scattering is based on the low-energy expansion 
of the 7rN Lagrangian in HBCHPT 

^ + £^ = 4 2 i + 4 4 7 ) + ... + fS + S + £S + ... . (i) 

The pion isotriplet field $ is represented by the field u(x) or U(x), where 
U = u 2 . In the so-called sigma parametrization 



where r are Pauli matrices and F is a LEC of the meson Lagrangian C^l (pion 
decay constant in the chiral limit). The nucleon isodoublet field \& = is 
decomposed using a time-like unit 4-vector v 

N v (x) = exp[imv ■ x]P+ty(x) 

H v (x) = exp[imv ■ x]P~ty(x) (3) 



v 
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(where m is the nucleon mass in the chiral limit) and the "heavy component" H v 



is integrated out |L2 



The Lagrangian ([!]) was constructed from the following building blocks: 

Up = 2-{V(<9 M - ir^u- u(d^- tt^v)} 
T M = - ir^u + u(d^ - ii^)u ] } 

X± = 2B{v)(s + ip)v) ±u(s + ipYu} (4) 



where B is a LEC of the meson Lagrangian 4 2 2 and s, p, £», r M , # are external 
fields (scalar and pseudoscalar, left-handed and right-handed vector isotriplet and 
vector isosinglet, respectively). For the elastic 7rN scattering one can set these 
fields to zero (dots correspond to terms vanishing in such a case) with the only 
exception of the scalar field. Via this field the nonzero quark masses are taken 
into account, and in the isospin limit[]m M = one has 

2Bs = M 2 (5) 

where M is the bare pion mass. 



1 From now on we shall work within the isospin limit, without stating it always explicitly. 
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Table 1: The field monomials Pi in £^ contributing to the elastic 7rN scattering. 
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Table 2: The field monomials Oi in C^ N contributing to the elastic 7rN scattering. 



Particular terms in ([!]) are |3[ [13| 



£(2) 
£(4) 



(i) 

-iV 
(2) 



E i(u-u + x+) 
— TIP 

N v (iv-V + g A S ■ u) N v 



N, 



ill 



(V • V + ig A {S • V, t> • «}) + ax(u ■ u) 



+ a 2 ((v ■ u) 2 ) + a 3 ( x+ ) + a 5 ie^ a v p S a UpU u + ...)#« 



(6) 



a 5 - 



i-3# 



+ |S^!i-V + ^{w • u, u^ie^VpSaiVu + h.c. 



where = ifi^o^v v is the spin matrix, g A is the neutron decay constant in the 
chiral limit, (.) denotes trace, and dots correspond to terms not contributing to 
the elastic 7rN scattering in the isospin limit. The field monomials and Oi 
together with their y+ and /3, (defined below) are collected in Tables [3] and ||. 
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Some of the LECs are divergent and are decomposed in the standard way 
(using dimensional regularization) as 



h = ^ r 0)+7i£0) m 
h = 6K/^) + (4vr) 2 AL(/i) {l) 

where 

m = g?{^i-|[ln47r + l+r'(l)]} . (8) 
The scale dependence of II (fJ,) and is given by 



FM = /[M + ^ln^f 



(9) 



(10) 



In what follows we shall often use /xq = M n and express results in terms of U and 
hi defined by 

Finally let us note that the nucleon field N v used in (|S|) is not exactly the 
one defined in (|3[), but is rather obtained from (|3D by so-called EOM (equation 
of motion) transformations fll3 |. However these EOM transformations have no 
effect on S-matrix elements. 



2.2 Renormalization 

To express the bare constants in C^} and in terms of measurable quantities, 
one has to calculate one loop corrections (using these lowest order Lagrangians) 
to the pion and nucleon propagators, as well as to the coupling of the external 
axial field a M = (r M — i^)/2 to a pion and to two nucleons. Tree contributions 
from the higher-order Lagrangians should be included as well. The results for 
pions are 

M 2 = Ml{l + ^h) 

= 1-f (2ll + 6m + ^\nf) (11) 

f = F^i-^q 

and for a nucleon with momentum mv + k 
m = m^v — 5m 

Mk) = i-S + I ^ £ -W( 6L ^ + ^ + ^ ln t) (12) 

9 a = 9A(l + ^) + T^(gl-4bu) 



where 



AMI gg^g run 
dm = a 3 - . 13 
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Let us comment on the results for nucleons, since they are different from what 
was used in most of previous calculations (see e.g. ||14j| ). The first difference is 
the term Sm/rriN in Z^- factor and cja- The reason of this difference is, as pointed 
out in ||15|| , that the Lagrangian commonly used in previous calculations contains 
some terms, which were transformed away in (|(J) by EOM transformations [ 13] . 
The second difference is the term [k — 5m.v) 2 j 'Am^ in Z/v-factor, which accounts 
for the contribution of heavy component of the nucleon source. This issue was 
not discussed in the literature on HBCHPT yet, so we shall briefly comment on 
it. More systematic analysis is to be found in [16| . 

In the path integral derivation of HBCHPT Lagrangian [Q , one starts from 
the generating functional of Green functions with at most two nucleons 



Z[r),rj\ 



exp i js™ + S nN + J d 4 x (rpf! + #77) J 



(14) 



After decomposition of the nucleon field to heavy and light component fl3j) and 
similar decomposition of the nucleon source 



one writes 



Pv(x) = exp[imv ■ x}P^~rj(x) 
Rv{%) = exp[imv ■ x}P~rj(x) 



S nN = J d A x [N V AN V + H V BN V + N V B'H V - H V CH V } 



(15) 



(16) 



where B' = ^B^ . After Gaussian integration over 



VH„VH, 



one gets 



Z[p v , R Vt p v , -R, 

SttN 



VNVNVu 



A 



SttN + sources 



/ d A x C nN = J d 4 x N v {A + B'C^B) N v 



;i7) 



sources 



/ d 4 x [(p v + RvC-'B) N v + N v ( Pv + B'C^Rv) + R v C~ 1 R l 



where the determinant Ah turns out to be just a constant. 

In S-matrix elements, the source rj is, so to say, replaced by a Dirac bispinor 
u(mv + k, a) (for a nucleon with momentum mv + k and spin a). The corre- 
sponding light and heavy components lead to 



Pv 

R v 



P+u(mv + k, a) 
P~u(mv + k, a) 



i 



f 1 4- 

V ~ 2m N 



2mjv 



2mi\ 



u(mv + k, a) 



u(mv + k, a) 



Clearly, p v leads to a quantity of the zeroth, and R v of the first chiral order. The 
lowest chiral order in the products C~ 1 B and B'C~ X is also the first one. So 
the heavy source terms in fllTl) contribute two chiral orders higher than the light 
source ones, and they start to contribute at the third chiral order. 
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To incorporate heavy sources into calculations of various amplitudes does not 
require any extra effort. For every Feynman diagram with light sources, one just 
adds diagrams with one or both light sources replaced by the heavy ones - the only 
differences in corresponding amplitudes are different factors for nucleon external 
legs. As a matter of fact, the whole effect of heavy sources is finally reduced to 
the momentum dependent, but otherwise very simple term in Z/v-factor. Using 
Z N given by fll2j) , one can calculate with just light sources and forget completely 



about heavy ones, since all their effect is taken into account in this JT6 . 

Heavy sources were not considered explicitly in most of previous calculations. 
The reason is that in the third order they only appear as the Zjy corrections to 
the first order, and the first order used to be treated in a special way — namely 
one used the original "relativistic" Lagrangian C^ N and then expanded the result 
in \i — M n /rriN. Advantage of this approach is that it gives automatically not 
only the lowest order amplitude, but also higher order corrections to it, including 
the discussed part of the Z N correction. Disadvantage, from our point of view, is 
that part of these higher order corrections is included in the definition of LECs of 
higher orders. So to avoid double counting in "relativistic 1-st order" approach, 
one has to know precisely what portion of LECs comes from the lowest order 
Lagrangian. This is well known for the second order LECs, but not for the third 
order ones, and if one does not want to work this out, it is better to perform the 
whole calculation strictly within HBCHPT. 

However, we can use "relativistic 1-st order" approach to check our prescrip- 
tion for the treatment of heavy sources. Let us take e.g. Born term for the elastic 
7rN-scattering (i.e. first order term with the nucleon propagator). "Relativistic 
1-st order" calculation gives for the invariant amplitude A + defined below in (|35|) 

At,= 9 -^f ■ (19) 

On the other hand, HBCHPT calculation of the same term in tree approximation 
up to the third order, ignoring heavy sources and expanding in /z, gives (for a 
kinematics defined by fl2"3|), (p?|)) 



A+ = ^(m N --^-)+... (20) 



Fl \ 8m N , 

and heavy source contribution, which turns out to give just an overall multiplica- 
tive factor (1 + t/&m 2 N ) , brings these two results in agreement. For the invariant 
amplitude B + one has the same situation, since also here one gets 

B+ = B ^( 1 -sk) + - ■ (21 > 

For the invariant amplitudes A~ and B~ , "relativistic 1-st order" and HBCHPT 
results agree to order 0(fi 2 ) even without heavy sources contribution, and this 
agreement is not spoiled by the heavy sources correction, because it starts to 
contribute only at order C(/i 3 ). 
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3 Elastic 7rN scattering 
3.1 Kinematics 

For the elastic 7rN scattering we use the following kinematics: ingoing pion with 
momentum q and isospin a, outgoing pion with momentum q' and isospin b, 
ingoing nucleon with momentum mv + p, and outgoing nucleon with momentum 
mv + p' 

n a (q) + N(mv + p) -> ir\q') + N{mv + p') . (22) 

The natural (and advantageous) choice for the 4-vector v is the 4-velocity of the 
ingoing (on-shell) nucleon 

m N v = mv + p . (23) 

In such a case one has p = 5m.v, where Sm is the sum of small quantities of 
the second and third orders, and so p itself is also a small quantity of the second 
order. 

The scattering amplitude T is given in the form 

T = T + 5 ah - T-ie abc T c 

(24) 

T± = u N (p', a 1 ) (a± + ^ie^q^VpSu) u N (p, a) . 

where 

un(p,ct) = P+u(mv+p,cr) = u{mv + p,a) 

u N (p',a') = P+u(mv + p' , a') = [l + - p* )] u(mv +p', a') . 

(25) 

The kinematics enables a significant simplification in the highest order (third 
one in our case) of the chiral expansion of the scattering amplitude . It is based 
on the fact that for the on-shell outgoing nucleon, v ■ p' is of the chiral order 2 
(although p' = p + q — q' itself is of the chiral order 1). 

m% = (mv+ P r =► V -P' = !I ^^~^ = 0{p2) ■ (26) 

So one can neglect v ■ p' (as well as v • p) in the highest-order amplitude. Or in 
other words, one can replace every v ■ q' by v ■ q in the highest-order. The error 
thus introduced is of even higher chiral order, beyond the level of precision of the 
calculation at hand. 

For sake of brevity, the following notation is used w = v ■ q, w' = v ■ q' 
and if needed, one can express the amplitude in the usual Mandelstam variables 
s = {rriNV + q) 2 , t = (q — q') 2 and u = {m^v — q') 2 using 

q-q' = Ml-\ 

v-Q = 2^(*-™ 2 n-M 2 ) (27) 
W = ^(s-m 2 N -M 2 + t) . 
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3.2 Scattering amplitude 

We shall present the amplitude order by order in the chiral expansion 

T = + 7(2) + 7(3) . (28) 

7(i) is given by the tree graphs (Fig. 1) containing vertices defined by 
(labelled as "1-st" in Figures 1-3), with the bare coupling constants and masses 
(F, qa-i M and m) replaced by their physical values (F n , g^, M n and m N ) in the 
final formulae. This replacement (which will take place in all orders 7( n )) just 
means that after writing F = F n — SF n , etc. and expanding into a power series 
in the deltas (SF n , etc.), we keep in the given order of the chiral expansion only 
the lowest order in deltas, higher orders in deltas are incorporated into higher 
orders of the chiral expansion. (The reason for this standard procedure is that it 
guarantees exactly the same form of Tm in tree calculation, one-loop calculation, 
etc., and an analogous statement is true for any higher 7L).) 

7(2) is given by the tree graphs with one vertex from C^ N (labelled as " 2-nd" 

in Figures 1-3) and all the other vertices from C^ N (Fig. 2). Again the bare 
coupling constants and masses are replaced by their physical values. 

Finally, T (3) is given by the sum T^ c + T^° p + T^. Here T$ c is given by the 

tree graphs with either one vertex from C^ N (labelled as "3-rd" in Figures 1-3) 
or two vertices from Cfh, and all the other vertices from C^ N (Fig. 3). The loop 
contribution T}°? p is given by one-loop graphs with vertices from Cj^l (Fig. 

4). As always, the bare constants are replaced by their physical values. The third 
term is just the contribution from this replacement in lower orders. 

Both T%? e and Tj°? p are infinite and renormalization scale dependent. It i 



is 

preferable to have the result in the explicitly finite and renormalization scale inde- 
pendent form, so we shall shift all the infinite and /z-dependent terms from T^° p 
to T^ cc . These shifted terms are polynomials in external momenta (otherwise 

they could not be canceled by the counterterms from C^jj ) and it is reasonable 
to shift from T}°? p to 7^ cc also other possible polynomial terms (in external mo- 
menta, including terms with negative powers of w), so that after adding also 
to T^ cc , one has a clear separation of the amplitude into the part TS? 1 , polynomial 
in external momenta plus poles in w, and the rest T™, which contains imaginary 
part, required by unitarity. 



8 



(a) (b) 

Fig.l Feynman diagrams for nN scattering contributing to the first 
chiral order. Double and single lines correspond to nucleons and pions, 
respectively. Crossed diagrams are not shown. 




(a) (b) (c) (d) 

Fig. 2 Diagrams contributing to the second chiral order. Black circle 
is the second order counterterm. Crossed diagrams are not shown. 
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(s) (t) (u) (v) 



Fig. 4 1-loop diagrams contributing to the third chiral order. Crossed 
diagrams are not shown. 
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We shall therefore present results at the third order in the form 



Tm=T^ + T^ . (29) 

Tj-g™ is made finite and scale independent by hand, T^ 1 by the cancellations of 
divergent and scale dependent terms from loops and counterterms. The coun- 
terterms in the third-order Lagrangian (Table |2[) were completely fixed by the 
general renormalization of the generating functional |17j], and so this cancellation 
provides us with a nontrivial check of the calculations. 

Calculation of T from (|1|) is in principle straightforward. Here we present 
just the final result, details are postponed to appendix [(J. For loop functions the 
following notation is used 

Jl™(w) = J$ ni (w) ± J$ ni (-w) 

Jr(w) = J (w)+4wL(fi)-^(l-2lnf) (30) 
ir(t) = / (t) + 2L(^)- T ^(l-21n. 



M 
I- 1 



and the loop integrals Jo, Jo, and K Q (which will appear later) are displayed in 
appendix ||. All the loop integrals in the amplitude are understood with the bare 
pion mass M replaced by the physical mass M n . 

Let us remark that two of the LECs from Table |2| — namely 64 and hyj - 
do not appear in the result. The reason for 64 is that the term containing this 
LEC is of the form wp ■ (q + q') which is of the fourth chiral order, since p itself 
is of the second order for the on-shell nucleon. The reason for 6 17 is, on the other 
hand, cancellation (in the third order) of different terms containing this LEC. All 
the other LECs from Table |2| are present in the result given belowQ: 

4) = 



(2) mjv 



Yi (aiq -q' + a 2 ww' - a 3 M%) 



+ 8m N F% 

+ 6ww' — 2q ■ q' 



(31) 



, pol _ g 2 A M* , W 2N . g\Ml 1 / , 2 

"(3) - 32ttF4 l 4 ^ 9 6IV1 it) + 12ttF4«,2 KQ Q W 



4r = ^J^\w)-{^{2Ml-t){Ml-2t)K ^t) 

+ ^(gV -w 2 ){Ml-w 2 )J™{w) 



2 For sake of completeness let us note that when taking into account also isospin breaking 
effects, three more LECs would appear, namely CJ4, &is and 620- 
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^(1) ~~ 2F% \w ^ w'J 
pol _ 1 

(3) 



o+ Pol 



o+ uni 
P(3) 



4tt 2 F4 



(&16 - &15) ™ - &19 



+ (g-g'-W) + ^(™W) 

- ^) (Wg • g' - M* W - M 2 g • q' - w 2 w' 2 ) 
(6x + 62) g • g' + + 2& 6 M 2 - 6 19 ^ {q-q'-w 



8m : 



^ 1 »i 

-'2 "1" Ml 



3(M 2 -g-g')} 



5^*; [\%w 2 - (7 + 11^) M 2 - (5 + 13g 2 A ) q • g'j 
Ti £^i ; (g • q' - w 2 ) (6M 2 - 5w 2 ) + (M 2 a 3 + ±t) 



1 



?,F: 



W 



3w 2 + g\±(q-qi-w 2 )(M 2 -u 
(Ml-\t)+2g 2 A (M 2 -lt)]ir(t) 



A. 

IFl 



m N l 



(1-1.) 

\w w' J 



16vrF4 ~ 12yrFi 



& (t - AMI) ^o(0, t) + §tM M * - ™ 2 )^ 
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4 Comparison with data 



4.1 Partial waves and threshold parameters 

To compare with experiment, it is convenient to express the results in terms of 
v4 ± and B ± , defined by the standard parametrization of the elastic 7rN scattering 
amplitude 



T nN = u(mv+p',a') 



A ± (s,t,u) + B ± (s,t,u) 



u(mv + p, a) 



(35) 



Starting from (^) one derives (after some algebra ) following relations between 



of 



f3 ± and A ± , B ± (see also Ref. @) 



A ± 
B ± 



-m N + 



4mjv 







± 



(36) 



Experimental information about elastic 7rN-scattering is usually given in terms 
of partial waves. For low energies, partial waves are characterized by a small 
number of threshold parameters such as scattering lengths, volumes, effective 
ranges, etc. These threshold parameters are not directly measurable, they can 



be extrapolated, however, from the experimental data above threshold [18] [IS] . 
Since the closer to threshold one is, the better CHPT is expected to work, these 
parameters seem to be the most suitable quantities to calculate and to compare 
with the (extrapolated) experimental values. 

Partial wave amplitudes for the scattering of a spinless particle on a spin \ 



particle are given by [18 



^ J\{(E + m iV )[A ± ( S , 9) + (VS - m w )S ± ( S , Oncost 
(E — m N )[A ± (s, 9) - (v/i + m N )B ± (s,9)]Pi ± i (cos 9)} dcos9 



(37) 

where E is the nucleon energy in CMS, and 9 is the scattering angle in CMS. 
The low-energy behavior of the partial wave amplitudes is characterized by the 
threshold parameters af ± , bf± defined by the expansion 



Re/, 



i±\ 



at ± k» + bt ± k^ 



+ . . . 



(31 



where k is the 3-momentum of a particle in CMS. 

Calculation of these threshold parameters from the scattering amplitude (|31|- 
|3~3| ) is in principle straightforward, results for the first four partial waves are given 
in appendix 0. 
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up to 0(p) 


up to £>0 2 ) 


up to 0(p 3 ) 


exp. 




-55 


-59 


-36 


-36 ±7 


a 2+ 


55 


59 


56 


64 ±3 


4+ 


180 


210 


280 


440 ± 140 


Os- 


-30 


-31 


31 


160 ± 120 


a 3+ 


-180 


-210 


-210 


-260 ± 20 


a 3- 


30 


34 


57 


100 ± 20 



Table 3: Comparison of the HBCHPT results for two D-wave and four F-wave 
threshold parameters up to the first, second and third order, with (extrapolated) 
experimental values. Units for D-wave are GeV~ 5 , for F-wave they are GeV~ 7 . 



4.2 Threshold parameters independent of LECs 0$, bj 

The scattering amplitude (|3~TH34]), and therefore also the sixteen threshold pa- 
rameters given in appendix |B|, contain four LECs of the second-order 7rN La- 
grangian (aj, 03, 05) and five linear combinations of LECs of the third-order 
7rN Lagrangian {b\ + b 2 , 63, b e , &i 6 — 615, 619). However, six out of these six- 
teen threshold parameters do not depend on any of these LECs. So we have at 
our disposal six parameter-free (where by parameters the LECs Oj and bj are 
understood) HBCHPT predictions. 

To compare them with experimental data, one needs some consistent set of 
extrapolated threshold parameters, for which we take Ref. [p0|1 . This comparison, 
which is shown in Table ||, is neither too impressive, nor discouraging.^] Two main 
lessons can be drawn from Table U 

• There is a clear improvement, when going from lower to higher orders of 
HBCHPT. 

• Contributions of different orders are comparable, i.e. the chiral expansion 
converges slowly. 

Our results for D- and F-wave threshold parameters are in agreement with 
independent calculation of these parameters in | I0| . In this work the Born term 



was treated in different manner, namely the "relativistic 1-st order" approach 
(discussed above) was used and results are presented to the leading (F-waves) and 
next-to-leading (D-waves) order in powers of fi = M^jm^. Our calculation, on 
the other hand, is fully based on the HBCHPT Lagrangian. For every threshold 
parameter, the two approaches give the same result to the order in powers of \i 



given in fllQ| , but there are differences in higher orders of /x, leading in some cases 
to nonnegligible (but not at all dramatic) differences in numerical values of the 
threshold parameters. 

3 One should be aware that experimental errors in pp[ are underestimated, since they refer 
" only to that part which can be estimated from deviations from the internal consistency of the 
method" . 
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-2.60 ±0.03 


a 2 


1.40 ±0.05 


a 3 


-1.00 ±0.06 


a 5 


3.30 ±0.05 


h + b 2 


2.4 ±0.3 


h 


-2.8 ±0.6 


h 


1.4 ±0.3 


bw - 615 


6.1 ±0.6 


big 


-2.4 ±0.4 



Table 4: The values of LECs from the x 2 -fit to the 10 7rN scattering thresh- 
old parameters, the nucleon cr-term, and the Goldberger-Treiman discrepancy. 
Uncertainities refer to the parabolic errors of the MINUIT routine. 



4.3 Remaining threshold parameters and LECs bj 

Let us now turn our attention to the remaining ten threshold parameters. As 
mentioned above, they contain nine LECs aj, bj and/or their linear combinations. 
Two of them are present also in HBCHPT results for two other quantities closely 
related to 7rN scattering. The first one is the pion-nucleon cr-term 

m{uu + dd) P) (39) 



a = (P 

2m N \ 

where P is the nucleon 4- momentum, m = (m u + m ( i)/2 and u, d are quark fields. 
The second one is the Goldberger-Treiman discrepancy 1 — g A m N / 1 F^g^N ~- Up to 
the third order of HBCHPT [| 

a = _!M 03 _ MM (40) 

m N 64:71 

^-^ttI 1 -^ 6 ") ■ (41) 

Values of the LECs can now be fixed by fitting the HBCHPT results to the 
(extrapolated) experimental values of the remaining ten threshold parameters, 
cr-term and g^N- Results of such a procedure, using values M n = 138MeV, 
m N = 939MeV, F w = 93MeV, g A = 1.26, a = 45 ± 8MeV, g„ N = 13.4 ± 0.1 
and threshold parameters from |20| , are presented in Table ^. (Let us recall that 
values of b%e — 615 and 619 depend on the choice of mesonic Lagrangian Hl3| .) 

One should stress that errors given in Table ^ are probably too optimistic, 
as a consequence of probably too optimistic errors of the used (extrapolated) 
experimental values of the threshold parameters. Another consequence of these 
probably underestimated errors is the relatively large value of the cr-term, a = 
59 ± 5MeV (corresponding to 03 = —1.00 ± 0.06). Enhancing error bars of the 
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up to (9(p) 


up to 0(p 2 ) 


up to 0(p 3 ) 


exp. 


units 


at 

u 





-0.13 


-0.07 ±0.09 


-0.07 ±0.01 


Gel/" 1 


bt 

U 





-24 


-13.9 ±3.0 


-16.9 ±2.5 


Gel/" 3 


u 


0.55 


0.55 


0.67 ±0.10 


0.66 ±0.01 


GeV- 1 


bo 


7.2 


8.2 


5.5 ±6.7 


5.1 ± 2.3 


GeV~ 3 




17.6 


51.2 


50.4 ± 1.1 


50.5 ±0.5 


GeV~ 3 


Oi_ 


-35.3 


-1.7 


-21.6 ± 1.8 


-21.6 ±0.5 


GeV~ 3 


ar+ 


-16.9 


-29.3 


-31.0 ±0.8 


-31.0 ±0.6 


GeV~ 3 


a7_ 


-17.1 


6.7 


-4.5 ± 1.0 


-4.4 ±0.4 


GeV~ 3 




18.6 


4.6 


31.2 ±0.3 


44 ±7 


GeV~ b 


a 2 _ 


-8.8 


-17.6 


-5.0 ±0.2 


2 ±3 


GeV~ 5 



Table 5: HBCHPT predictions for the 10 elastic 7rN scattering threshold pa- 
rameters up to the first, second and third order, compared with (extrapolated) 
experimental values. Theoretical uncertainties are shown only in the third order. 



threshold parameters leads to a value of cr-term closer to the experimental one. 
Another possibility, how to avoid a large cr-term, is to fix directly from the 
cr-term and not to treat it as free parameter in the fitting procedure. When doing 
so, it turns out that the third order LECs bt are not sensitive to the precise value 
of the cr-term. Variation of 03 from —0.8 to —1.0 (corresponding to variation of 
cr from 45 to 60 MeV) leaves the third order LECs almost untouched. 

On the other hand, the third order LECs are sensitive to the value of g n N- For 
g wN = 13.0 ±0.1 one obtains 61 + b 2 = 3.3 ± 0.3, 63 = -3.7 ±0.6, b 6 = 1.4 ±0.3, 
b w - 6 1B = 7.9 ± 0.6 and 6 M = -1.0 ± 0.4 . 

The values of the second order LECs Oj are in good agreement with their 
recent determination (ai = -2.59 ± 0.2, a 2 = 1.57 ± 0.1, a 3 = -0.92 ± 0.1, 
05 = 3.48 ±0.05 in our conventions). These values were obtained from a fit to the 
cr-term and eight other observables, four of them being some linear combinations 
of the threshold parameters used in our fit and another four being subthreshold 
parameters, not used by us. So the agreement is not a trivial one. As a matter of 



fact, the quoted results of [ ID | were obtained in a fit with enhanced error bars for 
threshold and subthreshold parameters (for reasons mentioned above). Original 
error bars lead to a x = -2.68 ± 0.03, a 2 = 1.47 ± 0.02, a 3 = -1.02 ± 0.06, in even 
better agreement with our results. 

Once we have pinned down the LECs, we can compare explicitly our results 



for the threshold parameters with the values given in |20|| . This is done, order by 
order, in Table ||. 

At first sight, there is remarkable agreement between theory and experiment. 
This statement is however almost empty, since we have chosen the values of LECs 
just to obtain such an agreement. The nontrivial content of this agreement is that 
it was achieved with reasonable, i.e. not unnaturally large, values of LECs. 
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LECs 


rest 


units 


U 


0.09 


0.03 


GeV' 1 


bo 


-3.3 


0.5 


GeV~ 3 


a t+ 


6.7 


-7.5 


GeV~ 3 


Oi_ 


-13.4 


-6.5 


GeV~ 3 


<h+ 


-6.7 


5.1 


GeV~ 3 




-6.7 


-4.5 


GeV~ 3 




1.9 


24.7 


GeV~ 5 


a 2 _ 


1.9 


10.8 


GeV~ 5 



Table 6: Comparison of contributions from counterterms with third order LECs 
and contributions from loops + other counterterms (denoted together as "rest"). 
Comparison is made in the third order at the renormalization scale fi = lGeV. For 
the remaining 8 threshold parameters, not displayed in the Table, contributions 
from counterterms with third-order LECs vanish. 

A criterion for the natural size of LECs is that corresponding counterterm 
contributions are not too large when compared to other contributions (loops and 
counterterms with fixed coefficients and/or lower order LECs) of the same order. 
Of course, such a comparison requires choosing some typical renormalization scale 
/i. Numerical results for fi = lGeV are presented in Table ^|. Contributions of the 
counterterms with the third order LECs are comparable, or even considerably 
lower than the rest of the contributions to the third order. 

Finally let us remark that both lessons learned from Table |3] are confirmed 
by Table [5| — the chiral expansion converges to the experimental values, but the 
convergence seems to be rather slow, in a sense that contributions to different 
orders are comparable. This fact seems to show that despite of the relative success 
in describing elastic 7rN scattering at threshold, the third order is definitely not 
the whole story. A complete 1-loop calculation, which will include the fourth order 
of the chiral expansion, is probably needed for sufficiently reliable description of 
this process. 
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5 Conclusions 



We have calculated the elastic 7rN scattering amplitude in the isospin limit in 
the framework of HBCHPT, up to the third order. Since the chiral expansion is 
supposed to work well near threshold, we have used the extrapolated threshold 
parameters, like scattering lengths, volumes, effective ranges, etc., to compare 
the results with data. 

The elastic 7rN scattering amplitude and therefore also the threshold param- 
eters contain nine low energy constants (besides F n and qa from the lowest order 
Lagrangians) . All these LECs were fixed from the available pion-nucleon data 
- the pion-nucleon cx-term, Goldberger-Treiman discrepancy and the threshold 
parameters. Values of the second order LECs are in a quite good agreement 
with their recent determination in [|K| . Values of the third order LECs were not 
determined directly from ttN data until now. 

The third-order calculation brought a clear improvement in the description of 
data, and this improvement was achieved with naturally small LECs. The results, 
however, suggest importance of higher-order corrections, since the contributions 
of the first three orders are frequently comparable and this will probably be the 
case also in the fourth order. 
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A Feynman rules 



Propagators and vertices used in the calculation of elastic 7rN-scattering are sum- 
marized. Most of them can be found, e.g., in ||, we present them here mainly 
for the sake of completeness and also because of slightly different notation. Rules 
are given in the isospin limit m u = m d . Nucleon momentum in final state p' is 
outgoing, all the other momenta are ingoing. Momenta and isospins of pions are 
denoted as (qi, a) , (q 2 , b) , (q 3 , c) , (g 4 , d) . 



Propagators 

pion propagator 



is ab (42) 



nucleon propagator 



k 2 — M 2 + is 

(43) 

V ■ p + IE 



7rN vertices from C^ N 



one pion 



two pions 



three pions 



- y S.g r a (44) 



v-(qi ~ q2) e abc r c (45) 



AF 2 



2F 3 
four pions 



9 a 

[S-(qi + q2) SabT c + S.(qt + q 3 ) 5 ac r b + S.(q 2 + q 3 ) 5 bc r a ] (46) 



V.(qi + q2 + <?3 + qi) (Sabred + SacSbd + badhc)- (47) 
-§i?4 [v.(q 2 ~ qi) 5 cd £ abf T f + v.{q A - q 3 ) 5 ab e cdf r f + 
v-{q3 ~ qi) S bd e acf T f + v.{q A - q 2 ) 5 ac e bdf r f + 
w-(94 - Qi) 5 bc e adf T f + v.(q 3 - q 2 ) 5 ad e bcf T f ] 
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7rN vertices from C^ N 



counterterm 



one pion 



two pions 



./V AM 2 N 
\ 2m m 3 , 



9 A 
2mF 



v.q S.(p + p') r a 



Ai 



-^F2 (ai gi-92 + a 2 u.gi v-g 2 + a 3 M 2 ) 5 ab + 

[(9i - 92) • (p + p') + il6a 5 e'* ,, ' xr gi A1 g2i/Up5 , ff ] e a 6 C r c 



7rN vertices from £ 



(3) 



one pion 



two pions 



8m 2 F 
M 2 



q 2 S ■ q + 2 S ■ p q ■ p + 2 S ■ p q ■ p'^j r a + 



(619 - 26i 7 ) 5-gr a 



^4 {[-i64(u • 92 9i • (p + p') + u • 9i 92 • (p + p'))+ 
(616 - b 15 )e^ plJ q lfl q 2l/ v p S a v-(qi- q 2 ) + 

• 9i 92 M + u • 92 q^)e^{p - p') v v p S^ 5 ab + 

[-(h + b 2 ) qi-q 2 v- (q x - q 2 ) - 
b 3 v ■ q ± v ■ q 2 v ■ (q ± - q 2 ) + 2b 6 M 2 v • (g x - g 2 ) + 



7T 2 F 2 



^(8a 5 + 3 9 A -\)e^°q llx q 2v {p + p ) p S a 



£abc7~c 



} 



47r-vertex from CS 2 l 



jk{ [{qi + q2?-M 2 }5 ab 5 cd + 
[(91 + 9s) 2 - M 2 ]5 ac <5 M + 
+ g 4 ) 2 - M 2 ]5 a A c } 



pion counterterm vertex from 

M 2 



2% 



F 2 



(q 2 - M 2 )l A ~ M% 



Jab 
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= - L I 

1, J 



B Loop integrals 

Definitions 

A 

HQ 2 ) 
QMQ 2 ) 
g» v h{Q 2 ) + Q»Quh{Q 2 ) 

J (id) 

9^vh (w) + v IJ ,v u J 3 (lu) 
K (u,Q 2 ) 

g llu K 2 + Q fl Q u K 3 + 

{v^Qv + Qu,Vv)K 3 + 

{Qxg^v + Qf*gv\ + Qv9\»)K A + 
(vxg^u + v^g uX + v u g X/Ji )K' 4 + 
(vxQ^Qu + v^Q u Qx + v v Q x Q^)K b + 

(Q\Vlj,Vu + QpVvV\ + QvV x V )l )K' 5 + 

QxQ^QuK 6 + v x v^v u K' 6 



i r d D k i 



i J (2tt) d k 2 -M 2 



1 r d u k 1 

i J (2vr)" (k 2 -M 2 )((k+Q) 2 ~M 2 ) 

1 r d D k kjj, 

i J (2vr)" (fc 2 -M 2 )((fc+Q) 2 -M 2 ) 

1 r d D fc *V fe " 
i j (2it) u (k 2 -M 2 )((k+Q) 2 ~M 2 ) 



(2^ 


(k 2 - 


-M 2 )(ui- 


■v-k) 


d D k 








W D 


(k 2 - 


-M 2 )(uj- 


■vk) 


d D k 




k/j,ku 






(k 2 - 


-M 2 )(u>- 


-vk) 


dPk 






1 



i. J 





(k 2 - 


-M 2 )((k+Q) 2 - 


-M 2 )(u- 


-vk) 


dPk 




k 








(k 2 - 


-M 2 )((k+Q) 2 - 


-M 2 ){u>- 


-vk) 


d D k 












(k 2 - 


-M 2 )((k+Q) 2 - 


-M 2 )(u>- 


-v-k) 



(55) 



1 r dPk kxk^ky 

i J J^rj 17 (k 2 -M 2 )((k+Q) 2 ~M 2 )(ui-vk) 
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Results^ 



A = 2M 2 L(fi) 



1 M 2 
32^2 n ^ 



(56) 



-2^) + T^(l-lnf -rln 



1+r 



1-r 



[Q 2 < o] 

-2L(//) + (l - In ^ - 2r arctan ±) [0 < Q 2 < AM 2 } 

r In i±r + mr) [Q 2 > 4M 2 ] 



-2L( /U ) 



where 



16vr 2 



l-r 



1 - 



4Af 2 
Q 2 



(57) 



-4cL(/i) + & (l - In + arccosh ^ 



-AcoL(fM) + ^(l -In 



arccos - 



[u; < — M] 
[cu 2 < M 2 ] 



-4a;L(//) + ^ (l - In ^) - ( arccosh £ _ m ) [ u > M ] 



(58) 



^ (o,g 2 



arctan 



2M 



167TA 



In 



2M-y/Q 2 



16WQ 2 



In 



-2M 



[w • Q = and Q 2 < 0] 
[v • Q = and < Q 2 < AM 2 } 
+ m ) [v • Q = and Q 2 > AM 2 } 



(59) 

4 Divergent integrals are calculated using dimensional regularization, with the terms in de- 
nominators of integrals understood with the ie prescription from appendix A (i.e. fe s _ M 3 , v . k l _ u 
etc. in definitions are just abbreviations for k -i_^ t i +ie , v . k _^j +ie etc). Results are given in the 
limit D — ► 4, i.e. in the final formulae the dimension of space-time D = 4 — 2e is expanded in 
powers of e, then one uses = — jjp- + and afterwards e is sent to zero. 
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hm = -\h{Q 2 ) 

I 2 m = |((M 2 -¥)/ (g 2 )-iA + T ^(M^-f)) (60) 

W) = K(i-f)W 2 )-^A-^(f -I)) 

Ji (cu) = ujJq (uj) + A 

Mu) = i((M 2 -c 2 )JoM-cA + ^(cM 2 -l^)) (61) 
J 3 (a;) = Ji (a;) - J 2 M 

i^g 2 ) = -i^ (o,g 2 ) 
^l(o,g 2 ) = -/ (g 2 ) 

^ 2 (o,g 2 ) = |(M 2 -ig 2 )K (o,g 2 ) + iJo(o) 

^3(0, g 2 ) = (|-^)K (o,g 2 ) + ^Jo(o) 

^;(o,g 2 ) = -h(Q") 

^;'(o,g 2 ) = -^ 2 (o,g 2 ) 

^4(0, g 2 ) = -{JoW + lMZR^Q^ + lK^Q^ + lQiK^Q 2 ) 

K'MQ 2 ) = I - M 2 ) J (g 2 ) + |A + ^ (£ - M 2 ) 

^5(0, g 2 ) = -i/ (g 2 ) + ^A-^(o,g 2 ) 
ir;(o,g 2 ) = -^ 4 (o,g 2 ) 

ir 6 (o,g 2 ) = -^J (o) + ^i(o,g 2 )-^ 4 (o,g 2 ) 
^(o,g 2 ) = -/ 2 (g 2 )-3^(o,g 2 ) 

(62) 

where all the Kj(0, Q 2 ) are given only for v • Q = 
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C Feynman diagrams for 7rN scattering 

Amplitudes corresponding to Feynman diagrams in Fig.l - 4 are presented. On- 
shell conditions 

a 2 = a' 2 = M 2 
v ■ q — v ■ q = 0(p z ) y ' 

are used to simplify expressions. Amplitudes are given in terms of loop integrals 
defined in Appendix [FJ and the following notation is used 

J±{w) = J n {w)±J n {-w) . (64) 



first order 

Fig. la 



Fig. lb 



P" 
of 

P~ 



■ 2 

9 A 
AF 2 

2 
9 A 

2F 2 

.2 
9 A 

AF 2 
.2 
9 A 

2F 2 



a 

P" 
of 

P~ 



\w w' 
\w w' 
\w w' 

(11 

\w w' 








(q ■ q' 



1 f w + w') 



AF 2 





WW 



WW 



(65) 



(66) 



second order 

Fie. 2a 



Fig. 2b 



a 

P" 
a 

P' 









P + 




a 




P~ 


2 

_ 9 A 


4ww' 


SmF 2 


. 2 

_ 9 A 


iV _ 


AmF 2 
■ 2 

_ 9 A 


w 2 — 


SmF 2 


. 2 

_ 9 A 


|V + 
V w 


AmF 2 








. 2 

AmF 2 



g — (q ■ q' — ww ') 



. 2 

9 .JL- 



2mF 2 



W 



,'2 



W 



12 



+ (q.q'-2M 2 )(^ + 



(67) 



(68) 
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Fig. 2c 



fl+ 3a /J 1_\ 

l J " AmF 2 \w 2 w' 2 J 

0- _ g> 2 / 1 I 1 \ 

r " AmF 2 \w 2 ^ w' 2 J 



Fig. 2d 

" + = ^2 (aiq.q 1 + a 2 ww' - a 3 M 2 ) 

[3+ = 

a" = 

/5~ = -^2«5 



third order 



Fig. 3a+3b 



Fig. 3d 



+ 



2 
3a 



J_ 

4m 2 F 2 w 



■ 1 

9 A 



1_ 

8m 2 F 2 w 
■ 2 

9 A 1 

8m, 2 F 2 





^(26 17 -6 19 )-M^ + g-g' 



^ (26 17 - 6 19 ) - Ml] (q.q> 
(q-q -Ml + w 2 )(q-q -Ml 



Fig. 3c 



a 



+ 



P + = 
a~ = 
P~ = 



■ 2 

g A 



a 



+ 



P + = 
a~ = 
P~ = 



4tt 2 F 4 
1 

4vr 2 F 4 





J 16 



015 + 



(6i + b 2 )q-q' + hw 2 + 2b 6 M 2 
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Fig. 3e+3f 



P 

Fig. 4a, 4b 



a+ = 
P + = 



a 



p- = 

Fig. 3g 



" ~~ 6^rF 4 «,2W 9 W y 

/9 + = ^ (fir + 16a 3 

' 2 "f2 



«- = f&M^+ 16a 3)w-- 5 

34 M 3 1 
32ttF 4 u> 2 



Fig. 4c+4d 



Fig. 4e 



Fig. 4f 





■ 4 

9 A 1 

sf 4 "?^ 

■ 4 


(j+H - 


2J 2 (0)) (W-™ 2 ) 


/9 + 


9 A 1 
4pT^ 


(j 2 (w) + 


_ _^ M 2 W ) 


a - 


.4 

9 A 1 

8F 4 w 2 


(j 2 (w) + 




0" 


.4 

3.4 1 
4F 4 u> 2 


(j+(w) - 


2J 2 (0)) 



« + = ^^(M 2 - W 2 )J + H(g-g'- W 2 ) 

/3+ = |i^((M 2 - W 2 )J -H-2^A 

.4 



Fig. 4g, 4h 



P~ = 3 S^(M*-w*)J+(w) 

a + = (w 2 Jq(w) + 3wJi(w)j 

P+ = 

or = j^pi (w 2 Jq(w) + 3wJi(w)^ 

p- = 
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Fig. 4i 



w 2 



« + = ^(j 2 + H-2J 2 (0))(g-g' 

0~ = 4^( J 2 + H-2J 2 (o)) 



Fig. 4k 



a+ = 

P+ = 

oT = --^wl 2 (t) 

(3~ = 



Fig. 41 



a + = ^[(3t-M 2 )t^ 1 (0,t) + (llt-3M 2 ) J fr 2 (0,t) 

+ (5i - M 2 )tK 3 (0, t) + 10^(0, t) + 2t 2 fT 6 (0, t) + 6 J 2 (0)] 

/3+ = 

a" = ||w[tK;(0,t)+tir5(0,t) + (|-M 2 )/ (t) + iA 

r = -§4 2tf 2 (o,t) 



Fig. 4m 

a+ = 
P+ = 

/T = 

Fig. 4n, 4o 



Fig. 4p+4r 



or 

p- 



= 



■ 2 



= 



Fig. 4s 
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Fig. 4t+4u 



Fig. 4v 



of 
a 



+ 
















external legs (wave function) renormalization 



or 



§|i (±6Z N (p) + \6Z N (j/) + 6Z n 



2F- 



-W 



(lSZ N (p) + l5Z N (p')+5Z 7T ) 



p- = 




5Z N (p) 


4M%a 3 
m 2 N 


5Z N (p') 


AMlaz . 
m N . 


5Z n 


Ml 
8n 2 F% 



32tt 2 F2 



1 + 48vr 2 L(/i) + 3 In 



l + 487r 2 L(/i) + 3 In 



7 4 + 48vr 2 L( / u) + In ^ 
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D Scattering lengths, volumes, 



Threshold parameters for the first four partial waves are given in explicit form, 
order by order. In the third order, the general formulae are displayed only if 
they are reasonably short, otherwise only the numerical results are given (in 
appropriate powers of GeV). 

S-wave 



1st order 







K = o 



6b 



niNM^ (87) 
SnF 2 (m N +M^) 

2m 2 N +3M 2 -4 : g 2 i M 7v (m N +M 7v ) 
32nF 2 M^m N (m N +M W ) 



2nd order 



„+ {ai+a2-a3)Ml 
"0 - nFl(m N +M„) 



ai(4m%-2m N M n +M 2 )+a2(4m%+2m N M^+5M 2 )+a 3 (2m N M n -M 2 ) 
'0 ~~ 4nF$m 2 N (m N +Mn) 
g\ {2m 2 N +3m N M 7r +3Ml) 



6+ = 
+ 



(88) 



a = 
3rd order 

"jg\m N Ml 



256n 2 F*(m N +M 7r ) 

g\M^{l5Am 2 N -lSm, N M 7I +^Ml)-g\Ml(l2Srn N +&AM^) 
m72-K 2 F^m N {m N +M- w ) 



M%a 3 m lf M ^[l+fei+fe2+fe3+2b 6 ] 

°0 - 2-KF 2 m N (m N +M^) H 16vr 3 F*m N {m N +M„) 



(89) 



-0.16 + 1.48a 3 + 5.15 (&i + b 2 ) + 6.18& 3 + 3.826 6 + 0.646 19 
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-wave 
1st order 



a 



1+ 24nF%M v 



9\ 



a 



12n F$Mn 
M n -2g A (m N +M 7: ) 



1+ 487rF2M 7r (m JV +M 7r ) 

2m N M^-3M 2 -4g 2 A m N (m N +M 7T ) 
96irF 2 M^rriN (m N +M-n ) 



2nd order 



+ — 16m N ai+16M w a2+g A (m N +3M 7T ) 

~~ 48TTF 2 m N (m N +M 7r ) 

+ -a 1 (l6m%+l2M%)+a2(l6m N Mn-l2M%)+l2a 3 M%+g 2 A m N (m N +3M 7T ) 
~~ 487TF 2 m 2 N (m N +M w ) 

4a 5 (m J v+M^)+g^(m j v+3M 7r ) 
~~ 487rF2 mjv (m JV +M T ) 

_ 4a 5 (m ]v +M^)+^m iv 
a !- ~~ 2A-wF 2 m N {m N +M 7T ) 

3rd order 

= -6.47 + 0.62 (b 16 - b 15 ) - 0.78b 19 

a+_ = -8.60 - 1.24 (b w - b 15 ) + 1.566i 9 
ar+ = 2.61 - 1.08 (bi + b 2 ) + O.68&19 

a -_ = -6.94 - 0.014a 3 - 1.10 (&i + 6 2 ) - O.OI863 - 0.0356 6 + O.6861 
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-wave 



1st order 



h 2+ 



G0nF 2 M 2 m N 

g A (2m N +5M 7 y) 
480irF 2 M 2 m 2 N 



«2+ = 





60ttF 2 M 2 



m N 



a 2 _ = 



g 2 A (4m 2 N -6m N M 7r -WM 2 )+5M 2 
mOirFZMlm^imN+MTt) 



2nd order 



+ g 2 A (m N +2M n ) 

%2 + ~ 120nF^M n m 2 N (m N +M^) 

+ 80a 1 m N M n -80a 2 M 2 -g 2 A ^8rn 2 N +21m N M 7r +15M 2 ^ 
l 2- 960nF 2 M^m^(m N +M^) 



g A (m N +2M^) 
2 + " 120ttF 2 M w m 2 N (m N +M n ) 



a 2 - 



20a 5 M^(m N +M^)+g 2 A (l2m 2 N +9rn N M 7T -5M 2 ) 
960nF 2 M n m 3 N (m N +M^) 



3rd order 



<4+ = 



193g 2 A m N 



<4_ 



*2+ 



<<2- 



U5200n 2 F^M n (m N +M„) 

= 25.0 + 0.18 (&i 6 - 6i 5 ) -O.22619 

[l+g 2 A (7-5n)]m N -5ng 2 A M 7T 
~ 144007T 3 FjMn (m N +Mn) 

= 11.5 + 0.306 (61 + 6 2 ) - 0.1936ig 
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-wave 



1st order 



*3+ 



9\ 

U0TrF 2 M%m 2 N 



Ct 3 _ = 



840-wF 2 M%m 2 N 



a 3+ = 



9\ 



U0irF%M3 m 2 N 



do = 



840n F 2 M%m 2 N 



2nd order 



+ g\(2m N +3Mir) 

3+ - 280wF 2 M 2 m%(m N +M 7v ) 



(4_ 



g 2 A (4m 2 N -15m N M w -14M 2 ) 
WmnFZMlmjjimN+Mn) 



l 3+ 



g 2 A (2m N +3M n ) 
' 2807TF 2 M 2 m- i N (m N +M^) 



g 2 A (4m jy — Mtt ) 



3rd order 



«3+ = 



73g A rrtjy 

752640vr 2 F*M3 ( mjv +M W ) 

g3 L (2190m 2 v -9457M2) 
225792007r 2 F4M|rrt JV (rrt ] v+M^) 



fl 3+ = 



m N [2+g 2 A (18-77T)]-7wg 2 A M 7I 
470400tt 3 F*M| (mjv ) 



23.3 
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